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O . Abstract 

o . 

CsJ ■ The recent results on anisotropic flow in ultrarelativistic nuclear colli- 

er , sions along with recent methodical developments and achievements in 

D | the understanding of the phenomena, are reviewed. The emphasis is 

given to the elliptic flow results. 
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1 Introduction 

■ In recent years, the subject of anisotropic flow in ultrarelativistic nuclear col- 

lisions attracts an increasing attention of heavy ion community. One of the 
main reasons for that is the sensitivity of anisotropic flow, and in particular 
O elliptic flow to the evolution of the system at the very early times p|. 

Anisotropic flow is defined as azimuthal asymmetry in particle distribution 
with respect to the reaction plane (the plane spanned by the beam direction 
and the impact parameter). It is called flow for it is a collective phenomena, 
but it does not necessarily imply hydrodynamic flow. It is convenient to char- 
acterize the magnitude of this asymmetry using Fourier decomposition of the 
azimuthal distributions. Then the first harmonic describes so-called directed 
flow, and the second harmonic corresponds to elliptic flow; non-zero higher 
^ ■ harmonics can be also present in the distribution. The corresponding Fourier 

coefficients, v n are used to quantify the effect 0. 

The two reasons for anisotropic flow are the original asymmetry in the 
configuration space (non-central collisions !) and rescatterings. In of 
elliptic flow the initial "ellipticity" of the overlap zone is usually characterized 
by the quantity e = ((y 2 — x 2 ))/((y 2 + x 2 ), assuming the reaction plane being 
a;z-plane. With the system expansion the spatial anisotropy decreases. This is 
the reason for high sensitivity of elliptic flow to the evolution of the system in 
the very early times 0, [|, of the order of 2-5 fm/c, independent of the model. 

Due to the lack of space we will not discuss in detail all recent develop- 
ments regarding directed flow. Briefly mention a couple. In || a very interest- 
ing qualitative prediction is given: it is shown that the radial flow (isotropic 
expansion in the transverse plane) and an incomplete baryon stopping should 



lead to a "wiggle" in the rapidity dependence of baryon directed flow; v±(y) 
should change it's sign three times with rapidity! This effect was also ob- 
served in a hydro model || . Once found it would be a strong evidence for the 
space-momentum correlations caused by radial flow. In a recent paper the 
important question of the role of the momentum conservation in directed flow 
measurements are discussed. Concrete recommendations for the analysis have 
been worked out. 

Recently, many new results regarding elliptic flow have been obtained in 
all directions: experimental measurements, improving analysis methods, and 
theoretical understanding of the underlying physics of the phenomena. I will 
try to mention the most important results in all three directions, but concen- 
trate mostly on the last two questions. A more complete picture of the recent 
experimental results can be found in ||. 



2 Improving the methods 

A significant progress in theoretical description of anisotropic flow demands 
the accuracy in measurements. Thus the corresponding methods are evolving 
in the directions of improving both, the statistical uncertainties, and under- 
standing systematics in the measurements. For the first, improving on the 
statistical errors, we mention extensive use of proper weights (the best would 
bewoc v n (y,pt)), which leads to the improvement of the reaction plane resolu- 
tion by 10%-20% and reduction of the statistical errors, and using the scalar- 
product approach ||. In the scalar-product method the flow is given by: 

v n (y,Pt) = — 7= = , (1) 

2 J(Q a n ■ Q b n ) 



where u = (cos(n0), sin(n0)) is a unit vector associated with a particle of a 
given rapidity and transverse momentum, and Q are flow vectors for the "full" 
event and subevents "a" and "b" : 

Qn = E U - ( 2 ) 

i 

where sum is over all particles in an (sub)event. This method (which is also 
easy to implement and analyze) in addition to the flow angle takes into account 
the flow effects on the magnitude of the flow vector and as a result gives smaller 
statistical errors. 

In is not possible to determine the reaction plane in the collision directly. 
Therefore, any measurement of the anisotropy in particle production with re- 
spect to the reaction plane is based on the measurements of particle azimuthal 



correlations among themselves. Those correlations to a different degree (de- 
pending on what exactly is analyzed) include the contribution from the cor- 
relations that are not related to the orientation of the reaction plane (e.g. 
resonance decays), and often called non-flow contribution. For a reliable in- 
terpretation of the results the non-flow contribution should be estimated or, 
better, measured. 
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Fig. 1: STAR Centrality dependence 

of the correlation between subevent flow an- 
gles. 
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Fig. 2: Simulations |27|. The results 
from 2- (triangles), 4- (stars) and 6-particle 
(crosses) cumulants. 



As an example, we discuss the estimate of non-flow contribution in STAR 
elliptic flow measurements [10|. An important observation for that is on the 
centrality dependence of the non-flow effects. The azimuthal correlation be- 
tween two particles can be written as 



(«n,l<, 2 ) = ( 



e «n</>i e -in0 2 



(3) 



where n is the harmonic, and the average is taken over all pairs of particles. 
The S n represents the contribution to the pair correlation from non-flow effects. 
Then, the correlation between two subevent flow angles is 



(cos(2(^ 2 a) - #H)> « M sub {v 2 2 +g)cx M suh v 2 2 + g, 



(4) 



where M su & is the multiplicity for a sub-event. Here, we have taken into account 
that the strength of the non-flow correlations scale in inverse proportion to the 
multiplicity: g oc 1/M su b- What is important is that the non-flow contribution 
to (cos(2(\l4 a ' ) — ^2 ))) * s approximately independent of centrality. The typical 



shape of (cos(2(^2 — ^2'))), see > f° r example, Fig. 1, is peaked at mid-central 



4 a) -* 2 fe) ) 

events due to the fact that for peripheral collisions, M su b is small, and for 
central events, v 2 is small. In the estimates [IIJ of the systematic errors, the 



authors set the quantity g = 0.05. The justification for this value was the 
observation of similar correlations for the first and higher harmonics (it has 
been investigated up to the sixth harmonic). One could expect the non-flow 
contribution to be of similar order of magnitude for all these harmonics, and 
model simulations support this conclusion. Given the value g = 0.05, one 
simply estimates the contribution from non-flow effects to the measurement of 
v 2 from the plot of (cos(2( 1 I / 2^ — ^2 ))) using Eq.(£|) (see circle-point error-bars 
in Fig.6). The relative contribution of non-flow effects is largest for very central 
and very peripheral bins (where, the reaction plane resolution is smallest!). 

Anisotropic flow is a genuine multiparticle phenomena (which justifies use 
of the term collective flow). It means that if one considers many-particle cor- 
relations instead of two-particle correlations, the relative contribution of non- 
flow effects (due to few particle clusters) would decrease. Considering many- 
particle correlations, one has to subtract the contribution from correlations 
in the lower-order multiplets and use cumulants instead of simple correlation 
functions |TT|. For example, correlating four particles, one gets 



(Mn,lMn,2<, 3 <,4> = V n + 2 ' 2 ■ vj n + 25 n . (5) 

In this expression, two factors of "2" in front of the middle term correspond to 
the two ways of pairing (1,3) (2,4) and (1,4) (2,3) and account for the possibility 
to have non-flow effects in the first pair and flow correlations in the second pair 
and vice versa. The factor "2" in front of the last term is due to the two ways 
of pairing. The pure four-particle non-flow correlation is omitted from this 
expression (see below on the possible magnitude of such a contribution). If 
one subtracts from the expression (|5]) twice the square of the expression (Q), 
one is left with only the flow contributions 

((Wn,l M n,2<,3<,4>) = («n,lMn,2<, 3 <,4> ~ 2(w„,l<,2> 2 = ~V* , (6) 

where the notation ((...)) is used for the cumulant. A very elegant way of 
calculating cumulants in flow analysis with the help of the generating function 
is proposed in [pT| . The simulations [pTfl , see Fig. 2, confirm that using 4- 



particle cumulants reliably removes non-flow contributions. It also shows that 
even in the presence of genuine 4-particle correlations (due to clusters decaying 
into 4-particles, such were introduced into simulations) those correlations are 
combinatorially suppressed compared to real flow correlations and there is no 
real need for use of higher order cumulants in flow analysis. 

The high precision results available in modern high statistics and large ac- 
ceptance experiments become sensitive to another effect usually neglected in 



flow analysis, namely, event-by-event flow fluctuations. The latter can have 
two different origins: "real" flow fluctuations - fluctuations at fixed impact 



parameter and fixed multiplicity (see, for example, flTJ, |TJ|), and impact pa- 
rameter variations among events from the same centrality bin in a case where 
flow does not fluctuate at fixed impact parameter. Note that these fluctua- 
tions affects any kind of analysis, including the "standard" one based on pair 
correlations. The reason is that any flow measurements are based on corre- 
lations between particles, which are sensitive only to certain moments of the 
distribution in i^. In the pair correlation approach with the reaction plane 
determined from the second harmonic, the correlations are proportional to v 2 . 
Averaging over events gives (v 2 ), which in general is not equal to (v) 2 . The 
4-particle cumulant method involves the difference between 4-particle correla- 
tions and (twice) the square of the 2-particle correlations. It is usually assumed 
that this difference comes from non-flow correlations. Note, however, that this 
difference ((f 4 ) — (v 2 ) 2 7^ 0) could be due to flow fluctuations. Let us con- 
sider an example where the distribution in v is flat from v — to v — f max - 
Then, a simple calculation would lead to the ratio of the flow values from the 
standard 2-particle correlation method and 4-particle cumulants as large as 
{v 2 y/ 2 /(2(v 2 } 2 - {v 4 )) 1 ' 4 = 5 1 / 4 « 1.5. 

3 The physics of elliptic flow 

Many important developments in this area: better understanding of the trans- 
verse momentum dependence of anisotropic flow in low pt region with the help 



of the "blast wave" model [[14], [L5|, attempts |TJ[] to describe t>2 (p*) in high 
p t region accounting for the parton energy loss - "jet quenching" , calculation 
in a parton cascade model |19|], analysis of the anisotropies in Color Glass 



Condensate [13], and detailed analysis of the elliptic flow in the hydro mod- 



els H, |H |5 



A very interesting development in this field is an attempt to calculate 
elliptic flow in Color Glass Condensate - classical field approach to describe 
ultrarelativistic nuclear collisions. One of the important consequences in this 
approach [jl3| is strong event-by-event fluctuations in v 2. As it has been already 



discussed, such fluctuations would manifest themselves in the difference of the 
flow results derived from 2- and 4-particle correlations. 



3.1 Transverse momentum dependence 



I would like to mention here what is usually referred to as "hydro inspired", 
"blast wave" or "expanding shell", models. Such models consider particle 
production from a thermal source in a form of an expanding shell with the 
radial expansion velocity having some azimuthal modulation. The case of 



directed flow was discussed in jlEfl. That model was used to fit E877 data 



and gave quite reasonable results [16|. The model was generalized for the 



elliptic flow case in IT3. Later it was used to fit STAR data and was further 



generalized for the case of the elliptic shape shell [fl5 1 . In this approach: 



, = / o 2 ^0 fe cos(20 fe )/ 2 ( at )ir 1 (A)(l + 2 g2 cos(20 & )) 
V2[Pt) / o 2 ^0 b J o (a t )^i(A)(l + 2 S2 cos(20 b )) ' U 

where Io, J 2 , and K\ are modified Bessel functions, and where a t ((fib) = 
(pt/Tf) sinh(p(0b)), and f3 t {4>b) = i m t/Tf) cosh(p(0b)). The assumptions of 
this model are boost-invariant longitudinal expansion and freeze-out at con- 
stant temperature Tf on a thin shell, which expands with a transverse rapidity 
exhibiting a second harmonic azimuthal modulation, p(<pb) — Po + Pa cos(20b). 
Here, <fib is the azimuthal angle (measured with respect to the reaction plane) 
of the boost of the source element on the freeze-out hyper-surface |L4|], and 
Po and p a are the mean transverse expansion rapidity (vq = tanh(p )) and 
the amplitude of its azimuthal variation, respectively. In Fig. 3, the fit to the 
minimum-bias data with s 2 = is shown as the dotted lines. The relatively 
poor fit led the authors to introduce a spatially anisotropic freeze-out hyper- 
surface, with one extra parameter, s 2 , describing the variation in the azimuthal 
density of the source elements, cx 2s 2 cos(20b). This additional parameter leads 
to a good description of the data, shown as the solid lines in Fig. 1. A posi- 
tive value of the s 2 parameter would mean that there are more source elements 
moving in the direction of the reaction plane. The model predicts a specific de- 
pendence of the elliptic flow on the particle mass. This mass-dependent effect 
is larger for lower temperatures (Tf) and larger transverse rapidities (po)- 

The behavior of ^(pt) at large transverse momenta is also very interesting. 
One of the possibilities is that the anisotropy at such transverse momenta is 
due to path length dependent nuclear modification of the parton fragmenta- 
tion function (jet quenching). High p t parton produced in the direction of long 
axis of the overlapping region exhibits more inelastic (in addition to elastic) 
collisions than that emitted along the short axis. It results in smaller prob- 
ability to fragment into high p t hadron. The effect depends on the density 
of the media and thus the observed anisotropy could serve as a measure on 



this very density JT8[ (and features of the energy loss itself). The transverse 
momenta, where v 2 (pt) saturates could also help in understanding the origin 
of the particles in the region of 2-5 GeV/c: do they acquire their transverse 
momentum due to multiple scattering or they come from a fragmentation of 
even higher p t parton? The preliminary STAR data [P5| , Fig. 4, support the 
idea of flow saturation at high p t . 
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Fig. 3: STAR f|. v 2 {p t ) with fits (Eq. @) 
including s 2 parameter (solid lines) and 
without (dashed lines). 
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Fig. 4: STAR g§. v 2 (p t ) for charged par- 
ticles. 



3.2 Hydro and low density limits 

The values of elliptic flow measured at RHIC are comparable to that in hy- 
drodynamic models. There was a clear disagreement at lower energies. As it 
has been mentioned in the introduction, for elliptic flow one needs rescatter- 
ings. Denser the matter and more rescattering means higher elliptic flow. The 
current understanding is that in limit of zero mean free path, the hydro limit, 
one gets the largest possible values of elliptic flow. 



Interesting that the flow values obtained in parton cascade calculations [19 



at different transport opacities could in principle significantly exceed the flow 
values from hydro calculations. It raises a question about validity of the as- 
sumption that the largest values of flow can be reached in hydro model. Two 
lines shown in Fig, 6 correspond to the results of this model for two values 
of transport opacity, corresponding to 35 (upper line) and 13 times higher 
gluon density as given by HIJING model. Note, however, that this opacity 
was calculating assuming 1 mb gluon transport cross section and assuming 
the hadronization picture when the number of gluons equals the number of 
hadrons. If one would, for example, consider a system as a constituent quark 
gas, one would have to increase the cross section approximately by 3 times and 



quark density by approximately 2 to 3 times. That would give the opacities 
very close to that needed to describe the data. 
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Fig. 5: Preliminary NA49 results 
flow centrality dependence. 



Fig. 6: STAR |2j. v 2 centrality depen- 
dence from 2- and 4-particle correlations. 
Lines are the predictions from a parton cas- 
cade model|l^] (see text). 
In the hydro limit elliptic flow is basically proportional to the original 
spatial ellipticity of the nuclear overlapping region [p], |23|, p4fl, v 2 oc e. In 
the opposite limit, usually called the low density limit [T2(], 21], elliptic flow 
depends also on the particle density in the transverse plane: v 2 oc e dN/dy /S, 
where 5* is the area of the overlapping zone. It results in a different centrality 
dependencies of the elliptic flow in these two limits. The comparison of the 

In this 



21 



results on elliptic flow from this point of view was first done in 
picture, the transition to the deconfinement would lead to some wiggles in v-ijs 
dependence, 
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Fig. 8: 1)2/5 as a function of particle ra- 
pidity density for three colliding energies. 
Hydro limits taken from fl23{ . 
One indication that at RHIC energies flow is still proportional to the par- 
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Fig. 7: PHOBOS results on V2(rj) overlayed 
with dN/dr) f26|. 



ticle density can be seen from Fig. 7 (taken from [26]) which show that 1*2(77) 
closely follows dN c h/dr]. The 3D hydro calculations [25| also cannot describe 
the pseudorapidity dependence of elliptic flow, once more indicating that the 
hydro description could be not correct in spite of the large values of t>2 mea- 
sured at RHIC. 

Taking the recent (year 1, ^snn = 130 GeV) STAR results on elliptic 
flow from 4-particle cumulants |27j (no non-flow effects) and plotting them 
along preliminary NA49 || and E877 results, Fig. 8, also suggests that even 
at RHIC energies elliptic flow continue to rise with particle density. (Note 
possible systematic errors in Fig. 8 of the order of 10-20% due to uncertain- 
ties in centrality measurements. However, this uncertainty does not alter the 
general trend). Decisive measurements could be already at full RHIC ener- 
gies of y/SNN = 200 GeV. If t>2 would continue to rise with particle density 
(which increases for about 15% in this energy range), it could give difficulties 
to hydrodynamic interpretation. 
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